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Following an extremely interesting idea published long ago, the work function at the outer 
crust region of a strongly magnetized neutron star is obtained using relativistic version of Thomas- 
. . . Fermi type model. In the present scenario, the work function becomes anisotropic; the longitudinal 

0^ ' part is an increasing function of magnetic field strength, whereas the transverse part diverges. An 

' approximate estimate of the electron density in the magnetosphere due to field emission and photo 

, emission current, from the polar cap region are obtained. 
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, The study of the formation of plasma in a pulsar magnetosphere is a quite old but still an unresolved astrophysical 
' issue. In the formation of magnetosphere plasma, it is generally assumed that there must be an initial high energy 
, electron flux from the magnetized neutron star. At the poles of a neutron star the emitted charged particles flow 
only along the magnetic field lines. Further a rotating magnetized neutron star generates extremely high electro- 
I static potential difference near the poles. This potential difference is the driving force and plays the major role in 
HH ■ the extraction of electrons as field emission or what is also called cold emission, from the crustal matter of strongly 
^"i _ magnetized neutron stars near the poles. The flow of high energy electrons along the direction of magnetic field lines 
. and their penetration through the light cylinder is pictured with the current carrying conductors. Naturally, if the 
conductor is broken near the pulsar surface the entire potential difference will be developed across a thin gap, called 
Q , polar gap. This is based on the assumption that above a critical height, from the polar cap, because of high electrical 
%^ conductivity of the plasma, the electric field (= Eq in this article), parallel with the magnetic field near the poles 
c/3 , is quenched. Further, a steady acceleration of electrons originating at the polar region of neutron stars, travelling 
along the field lines, will produce magnetically convertible curvature 7-rays. If these curvature 7-ray photons have 
[ energies > 2meC^, then pairs of e~ — will be produced in enormous amount with very high efficiency near the polar 
^-H ■ gap. These produced e~ — e"*" pairs form what is known as the magnetospheric plasma 0, S i5|, i6|, i7|, @| . 
J> The process of extracting electrons from the outer crust region of strongly magnetized neutron stars, including the 

most exotic stellar objects, the magnetars, one requires a more or less exact description of the structure of matter in 
^ that region j^, ,10,, JTJ . From the knowledge of structural deformation of atoms in strong magnetic field; the departure 
from spherical nature to a cigar shape, allows us to assume that the atoms in the outer crust region, which are fully 
ionized because of high density, may be replaced by Wigner-Seitz type cells of approximately cylindrical in structure. 
We further assume that the electron gas inside the cells are strongly degenerate and are at zero temperature. It 
is well known that the presence of extraordinarily large magnetic field not only distorts the crystalline structure of 
dense metallic iron, also significantly modifies the electrical properties of such matter. As for example, the electrical 
conductivity, which is otherwise isotropic, becomes highly anisotropic in presence of strong quantizing magnetic field. 
In presence of strong magnetic field iron crystal is highly conducting in the direction parallel to the magnetic field, 
whereas flow of current in the perpendicular direction is severely inhibited. The aim of this letter is to show that the 
5J] ' work function, which is the most important parameter associated with the emission of electrons from the polar region 
' of strongly magnetized neutron stars, will also show anisotropy in presence of strong magnetic field. In this article 
we have shown that the work function, associated with the emission of electrons along the field lines increases with 
magnetic field strength. Whereas its transverse component diverges, irrespective of the dimension of the cylindrically 
deformed atoms. The scenario is very much analogous to the charge transport mechanism in presence of strong 
magnetic field. To the best of our knowledge, the study of anisotropic nature of work function in presence of strong 
quantizing magnetic field, which has relevance, specially in the case of magnetized neutron star crustal region has not 
been studied before. It is also well known that the most important surface emission processes are thermal emission, 
may be enhanced by Schottky effect, field emission, caused by strong electric field at the poles and perhaps the other 
important process is the photoelectric emission induced by high energy cosmic photons. The last process of electron 
emission, if at all possible, has not been studied in this context. 

In this letter, following Q we shall obtain work function associated with the emission of electrons, along and 
transverse to the direction of magnetic field lines, and show that the later component is infinitely large in this model. 
In [l| the author have obtained work function related to emission of electrons from cold cathodes. It was associated 
with some engineering problems. We assume cylindrically deformed atoms in the crustal region with their curved 
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surfaces parallel to the neutron star surface. Which further means, that locally, the field lines and axes of the cylinders 
are parallel with each other. To obtain the work function for both the cases, let us consider fig.(l). Here the magnetic 
field is along z-axis and r-axis is orthogonal to the direction of magnetic field lines. To get the longitudinal part of 
the work function, we assume that an electron has come out from the cylindrically distorted atom, through one of 
the plane surfaces, and is at the point P. The electro-static potential at Q produced by this electron is given by 

where AP — PD — r and CQ = z. This equation can also be expressed in the fohowing form 

for z > or < h (2) 

For the charge distribution within the distorted atom, the Poisson's equation is given by V^(/) = Anen^ where ip is the 
electro-static potential and is the electron density, given by 

"^ = |§ E(2-'5.o)(M'-m'-2^/eB)l/2 (3) 

where v is the Landau quantum number with I'max the upper limit. The first factor within the sum takes care of 
singly degenerate 0th Landau level and doubly degenerate other levels with i/ 7^ and fj, is the chemical potential for 
the electrons, given by = Ep — e(j){r) = constant, this is the so called Thomas-Fermi condition, Ep is the Fermi 
energy for the electrons (throughout this article we have taken h = c = 1). Hence we can re- write the Poisson's 
equation in the following form 

= y (2 - 5,o)[(m + ecl>f -ml- 2veB\^l^ (4) 

l/=0 

To get an analytical solution for in cylindrical coordinate system (r, 0, z), with azimuthal symmetry, we put, nie = 0, 
i.e., the kinetic energy is assumed to be high enough and Vmax = 0, which is a valid approximation, provided the 
magnetic field is too high, so that the electrons occupy only their 0th Landau level. Now defining tp — /j, + ecj), the 
solution in the cylindrical coordinate system is given by 

/>oo 

with z < (5) 

Here = 2e^B/Tr and a(^) is some unknown spectral function. Now following [H, we assume that there exist a 
fictitious secondary field in vacuum, the work function is defined as the work done by this field in pulling out an 
electron from the material surface to infinity. The secondary field is expressed in coherent with ^'■^^(r, z) and 4){r, z), 
and is given by 

POO 

</.(^)(r,z) = e / M^r)f{0 exp(-e^K 

with z > (6) 

where /(^) is again some unknown spectral function. Then following using the continuity conditions of tangential 
and transverse components of electric field and the displacement vector respectively, we finally get 

Wf ^ -e' ^ - ( ^^X^' m^e' (7) 

where iJc*^"* « 4.43 x 10^'^G, the critical field strength for electrons to populate their Landau levels in the relativistic 
scenario. This equation gives the variation of work function with the strength of magnetic field (~ B^^^) associated 
with the emission of electrons along the direction of magnetic field. From this expression it is also obvious that 
for a given magnetic field strength, if me is replaced by uip, the proton mass or m/, mass of the ions, then, since 
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mi ^ vrip ^ me, the same kind of inequalities are also valid for the work functions associated with the emissions 
along ^-direction. Then obviously, it needs high temperature for thermo-ionic emission of protons or ions, high electric 
field for their field emissions and high frequency incident photons are essential for the photo-emission of these heavier 
components. 

Let us now consider the emission of electrons in the transverse direction. As shown in fig.(l), p is the position of 
an electron, came out through the curved surface, then the electro-static potential at q is given by 



where ap = ro, cq = r and pd = z. It can also be expressed as 



ro)2]i/2 



(8) 



for 2 > or z < 



(9) 



whereas (f>{r, z) and (r, z) remain unchanged. Then as before, from the continuity conditions on the curved surface, 
and using some elementary formulas for the Bessel functions, we finally get 
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where F = Ji[(i?-ro)C]/Ji(i?C) and G = Ja[{R-ro)(]/ Jq{R^). Putting /(^) in </.(^)(r, z) and following [l|, evaluating 
the trivial ro-integral and finally using some of the properties of Bessel functions, we have the work function in the 
case of emission of electrons in the transverse direction 
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This is a diverging integral for both R — > oo, or R — > 0, i.e., the value of the work function is independent of the 
transverse dimension of the cylindrical cell. For R — > 0, the second term on the numerator vanishes, whereas the 
rest can be integrated analytically and found to be diverging in the upper limit. On the other hand, for R — > oo, 
the second term becomes ^ tan(^i? — tt/A) (asymptotically). Which itself also diverges for certain fixed values of the 
argument, but dose not cancel with the infinity from the other part of this integral. Therefore in this limit also the 
work function becomes infinity. Although we have considered here the extreme case of ultra-strong magnetic field, we 
do expect that even for low and moderate field values, the work function corresponding to the electron emission in the 
transverse direction will be several orders of magnitudes larger than the corresponding longitudinal values. In fig. (2) 
the monotonically increasing straight line curve indicates the variation of Wf with B. This graph clearly shows that 
for low field values, Wf is a few ev in magnitude, which is of the same orders of magnitudes with the experimentally 
known values for laboratory metals. However, in this model calculation, since the magnetic field is assumed to be 
extremely high, we are unable to extend our calculation for very low field values. Also, in this model we can not show 
the variation of Wf from one kind of metal to another. 

As an application of this model we now obtain the density of ultra-relativistic electrons near the poles populated 
through field emission process. It was first introduced by Fowler and Nordheim [T^ to explain cold cathode emission 
phenomenon. The underline process is the tunneling of metallic electrons through surface potential barrier. Following 
Fowler and Nordheim we have the electron number density 
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where a is the fine structure constant, Ep is the electron Fermi energy and Eq is the field at the poles causing electron 
field emission. In the numerical evaluation of field electron density, we have taken both a constant E^^ (lO^^V/m), as 
given in the literature [2] and also as a function of magnetic field strength at the poles, given in In fig. (3) we have 
shown the variation of field electron density with the magnetic field strength. Monotonically increasing curves are for 
magnetic field dependent Eq. The upper, middle and the lower curves are respectively for Ue = 10~^, 10~^ and 10"'^ 



4 



times normal nuclear density. Whereas the decreasing curves are for constant Eq, with the same kind of variation for 
electron density. 

Next, we calculate the photo-electric current in the longitudinal direction. In the relativistic case, in presence of 
strong quantizing magnetic field at zero temperature the photo-electric current J is given by 

eB 

J = ^{Ef + ~ Wf)[l + l^maxit^max + 1)] (13) 

where E-y is the incident cosmic ray photon energy. In fig. (2) we have shown the variation of photo-electric current 
density with the strength of magnetic field. For the sake of illustration, we have taken Ej = IGeV. The variation 
for He = 10"'^, lO"** and 10~^ times normal nuclear density are shown by upper, middle and the lower monotonically 
decreasing curves respectively. If the photon energy is increased, then it is obvious that J will increase linearly with 
Ej for constant B and (or Ep). 

In conclusion, we mention that the main purpose of this letter is to show that (i) in presence of strong quantizing 
magnetic field the work function becomes anisotropic, (ii) the transverse part is infinitely large, (iii) the longitudinal 
part increases with the strength of magnetic field, (iv) low field values of work function are more or less consistent with 
the tabulated values and finally, (v) the use of longitudinal part of work function gives consistent result associated 
with magneto sphere electron density, reported in the published literature. In the present work, to the best of our 
knowledge, the anisotropic nature of work function in presence of strong quantizing magnetic field is predicted for the 
first time. Further, the diverging character of work function associated with the electron emission in the transverse 
direction is also obtained for the first time, and so far our knowledge is concerned, it has not been reported in the 
published work. We have also noticed that in the low magnetic field limit (within the limitation of this model) the 
numerical values of work function are of the same order of magnitude with the known laboratory data. However, we 
are not able to compare our result with the variation from one metal to another. As an application the use of magnetic 
field dependent work function to field emission gives results which are consistent with the published theoretical results. 

Since the results presented in this article are very much approximate in nature, we therefore conclude with a 
beqautiful quotation from a very old but extremely interesting article by Wigner and Bardeen [l^]- "It is perhaps not 
quite superfluous to have, in addition to a more exact calculations of a physical quantity, an approximate treatment 
which merely shows how the quantity in question is determined, and the lines along which a more exact calculation 
could be carried out. Such a treatment often leads to a simple formula by means of which the magnitude of the 
quantity may be readily determined" . 
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FIG. 1: Schamatic diagram for electron emission along and orthogonal to the direction of magnetic field. 
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FIG. 2: Variation of work function (Wf) with the magnetic field strength shown by monotonically increasing straight line curve. 
Variation of photo-electric current with magnetic field. Upper, middle and lower curves are respectively for 10~^, 10"'' and 
10~^ times normal nuclear density. 




FIG. 3: Variation of field electron density with magnetic field. Monotonically increasing curves are for magnetic field dependent 
electric field, whereas the set of decreasing curves are for the constant value of electric field as discussed in the text. For both 
the cases lower, middle and upper curves are respectively for 10~^, 10"'' and 10~® times normal nuclear density. 



